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I. INTRODUCTION 

Parametric amplification together with stimulated emission are fundamental 
mechanisms for generation of coherent radiation. Laser action is based on the 
latter whereas the former is the basis of optical parametric oscillations (OPOs). 
Even though, laser systems are encountered worldwide thanks to their intense 
commercialization, OPOs are, nowadays, more and more spreading (and not 
only in laboratories) owing, essentially, to their tunability on a wide range of 
frequencies. In fact, for a laser the emitted frequency is fixed once and for all to 
the manufacturing. In consequence, a pile of several lasers of diff"erent frequen- 
cies is needed when the operating system involves a whole range of frequencies. 
We thus measure the considerable advancement that OPOs bring to modern 
optical technology where large frequency variations are desired. Although the 
basic principles of optical parametric amplifiers (OPAs) have been known for 
over 40 years, OPOs developed quickly only in the last decade mostly for tech- 
nological reasons [1]. The use of periodically varying media greatly enriches 
the diversity of observable phenomena either in OPAs or in OPOs, allowing for 
existence of coherent field patterns, the study of which is the main goal of the 
present paper. 

Compared to lasers, OPOs have received much less attention in spite of their strong 

n 

interest both on the fundamental and on the technological sides [2|]. We recall that these are 
very frequency agile coherent sources with a wide range of possible applications including 
range finding, pollution monitoring and tunable frequency generation. They are also the 
key element for the production of twin photons and the realization of fundamental quantum 
optics experiments 3|. 

Basically an optical parametric amplifier (CPA) generates light via a three-wave mixing 
process in which a nonlinear crystal subjected to a strong radiation at frequency Up (pump 
beam) radiates two coherent fields at frequencies Ug (signal beam) and Ui (idler beam) such 
that the energy conservation law ujp = Ug+uji is satisfied. This energy conservation criterion 
may be interpreted in terms of photons where one photon at frequency Up is converted into 
two photons at frequencies Ug and Ui. This process is most efficient when the phase matching 
condition is fulfilled. It states that optical parametric amplification is favored when the three 
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interacting waves keep constant relative phases along their propagation inside the crystal. 
This implies kp = kg + ki, or equivalently np{up)ujp = Ug {ujg) ujg + rii (tUj) Ui where kj and rij 
are the wavevector and the refractive index at frequency Uj respectively (j = p, s, i). 

OPOs have recently appeared as physical systems whose modeling generates specific 
problems sharing common grounds with very general open questions such as those related 
to the appearance of complexity in spatially extended systems j4|. In fact, OPOs have 
become one of the most active fields in nonlinear optics not only for the richness in nonlinear 
dynamical behaviors [5] but also for thepotential applications of OPO devices including 
low-noise measurements and defection jfi, 7]. When driven by an intense external field, both 
OPAs; i.e. corresponding to a single pass through the quadratic crystal, and OPOs show 
a number of remarkable features among them are the localized structures (LS) or spatial 
solitons [8|, |9|]. These nonlinear solutions are generated in the hysteresis loop involving stable 
homogeneous states and the nonhomogeneous periodical branches of solutions. The latter 
that initiate the LS formation result from the modulational instability (often called Turing 
instability) and have been intensively studied in nonlinear optical cavities [l^. Stability of 
Ss and effects of spatial inhomogeneities on the LS dynamics have been addressed in Refs. 
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12|. 



Very recently another type of LS in the form of dissipative ring-shaped solitary waves 
have been generated in the regime where the steady state solutions are stable with respect to 
the modulational instability. Indeed, it has been shown in [l^ that OPOs can continuously 
generate spatially periodic dissipative solitons with an intrinsic wavelength. These mod- 
ulations spontaneously develop from localized perturbations of the unstable homogeneous 
steady state that separates the two stable states of an hysteresis cycle. This constitutes the 
counterpart of Turing spontaneous modulations initiated by extended perturbations. They 
occur in the wings of the 2D traveling flat top solitons (fronts, domain walls or kink-antikink) 
and give eventually rise to ring-shaped propagating dissipative solitons. Such non- Turing 
periodic dissipative solitons have also been predicted in liquid crystal light valve nonlinear 
optical cavity jl^ . 

In the more general context of optics the soliton ability to self-confine light bearp power 
makes it a promising candidate for future applications in information technology jl5| and 



image processing 



16 



171]. In this context, an issue of important significance is the devel- 



opment of a research activity aiming at taking advantage of more sophisticated geometries 



of cavities |3] and periodic modulations of the refractive index [ol, Q to generate, control, 
and stabilize optical LS. The existence of quadratic gap solitons in periodic structures, hav- 
ing spatial extension much bigger than the period of the structure, and thus well described 
within the framework of the multiple scale expansion, was studied in ^. 



stabilization of LS are possible since periodic transverse modulations are 



The control and 
mown, with the 
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advent of photonic crystals, to profoundly affect spatial soliton properties 

In this paper we address the problem of the formation and the dynamics of localized 
structures in presence of transverse periodic modulations of the refractive index in both 
OPAs and OPOs, concentrating on field patterns whose spatial scales are of order of the 
period of the medium. This is actually the main distinction of our statement from the 
previous studies exploring either slow envelope LSs [9|] or the tight-binding limit leading to 
nonlinear lattice equations [2c3]. Indeed, optical parametric wave-mixing systems are one of 
optical devices that can operate either as a conservative (OPAs) or a dissipative (OPOs) 
system. We have analytically investigated the existence and stability of conservative LS (or 
solitons) that appear in OPAs and the dissipative periodic patterns occurring in a parametric 
oscillation regime (OPOs). We have first identified the simultaneous conditions required on 
index modulation together with the incident pump for stationary LS to exist by means of 
Bloch waves approach. It results from our analysis that in a degenerate configuration (the 
signal and idler are identical) signal and pump fields can be found in form of phase locked 
LS (real solutions). They are strongly stable and located at the minimum of the index 
periodic modulation. To extend the analytical analysis to higher nonlinear regimes we have 
performed numerical simulations that allowed characterizing the different types of stable LS 
that can be supported by the system under transverse periodic modulations. 

The paper is organized as follows. In Sec. II we briefly recall the two-wave model 
governing the spatio-temporal evolution of the slowly varying envelopes of the pump and the 
signal waves including diffraction and transverse periodic modulation of the refractive index. 
The general conditions required for the existence and stability of conservative LS appearing 
in OPAs, by means of Bloch waves approach, are given. Extension of Bloch waves approach 
to take into account dissipative effects stemming from cavity losses is carried out in Section 
III, for a degenerate OPO with transverse periodic modulations, where for construction of 
periodic solutions we use a kind of inverse engineering allowing us to determine the external 
field from a given structure of the nonlinear Bloch state. Concluding remarks are summarized 
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in the last section. 



II. LOCALIZED MODES IN PERIODICALLY MODULATED OPTICAL 

PARAMETRIC AMPLIFIERS 

A. The model 

We start with the model describing optical parametric amplification where we take into 
account diffraction and the transverse periodic change of the refraction index (n = n^°^ + 
j = 0,1). In the degenerate configuration, the signal and idler fields are identical 
leading to the resonant frequency conversion Uq = 2uJi where ujq and uoi are the pump 
and signal frequencies respectively. Considering beam propagation along the z-direction, 



the non 

system 



inear interaction of the two waves propagating in the crystal is governed by the 



2l|: 



dzao = -V^ao + ^^o(^)ao + iaje'^^' , (la) 
d^ai = iV'^ai + iei{f)ai + 2iaoaie~'^^'' , (lb) 

where ao, ai are the envelopes of the pump and the idler, respectively, r = {x,y), V stands 
for the transverse gradient, and ij = ujj[n^p]'^ /c^kj. Hereafter an overbar stands for the 
complex conjugation. 

We will be interested in periodically varying refractive indexes. Let us further specify the 
model considering Uj^^ to be independent on frequency, i.e. assuming 

e,{r} = 2e,{r}^e{r}. (2) 

This assumption is not essential for construction of the field patterns, reported below, and 
is introduced only for the sake of concreteness, as the approaches we will consider involve 
numerical simulations. Thus, the theory developed here is straightforwardly generalized 
to any relation between eo{r) and ei(r), provided that their spatial periods coincide (as it 
happens in a typical experimental situation). 

For perfect phase-matching, Ak = 0, the system ([1]), ([2]) can be written in the Hamiltonian 
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form dz<yj = iSH/Saj with the Hamiltonian 

1 



H 



iVaoT - |Vai| 



+ e \ao\ + -\ai 



(3) 



being an integral of "motion": dH/dz = 0. Another integral is given by the total power 
P = 2Pq + Pi, where we have introduced a notation Pj = J \aj\'^d'^f for the power of the 
j-th component, i.e. dP/dz = 0. In what follows, existence domain and properties of 
non-uniform stationary solutions will be investigated. 



B. General properties of stationary solutions 

For the sake of simplicity, we address in the present work the situation of y-independent 
spatial patterns, i.e. we assume n^p = n'p{x). Without loss of generality we can impose that 
rij'\x) is a vr-periodic function (this can be always achieved by proper renormalization of the 
coordinate units). Then it is convenient to introduce a periodic function V{x) = — 2£(x), 
V{x + tt) = V{x), and to define two linear operators 

and the respective eigenvalue problems 

^-M^li^) = ^^}{^- (5) 

Here i^u}{x) are the linear Bloch waves, with the index v = 0,1,2,... standing for the 
number of the allowed band and q designating the wavevector in the first Brillouin zone 
(BZ), g G [—1, 1]. In this paper we will also use the notations £^"|j_ for the lower ("— ") and 
the upper ("+") edges of the z/-th band, and A^/^ = £[,+i__ — £if+ for a u-th finite gap {Aq^ 
designating the semi-infinite gaps). 

For a general form of V{x), the introduced linear eigenvalue problems (jS]) represent Hill 



equations 



22l | and, as it becomes clear below, play crucial role in the systematic construction 



of the nonlinear localized modes (for a review see e.g. [231]). 

Let us first investigate the existence of the different types of solutions of the conservative 
system ([1]), which are homogeneous in z-direction. Hence, we seek for non-uniform {x- 
dependent) stationary solutions of the system in the form: 

ao = ao{x)e-'^^^'+''y^\ ai = al(x)e-'('^^^+'^-^)/^ (6) 
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where we have taken into account the matching conditions = k^o = 2kzi, ky = kyo = 2kyi 
and Ak = . Substituting the ansatz in Eqs. ([1]) we obtain 

nao = + V{x)ao - 2a?, (7a) 

—ai = --^ + ^[x)ai-2aQai, (7b) 

where VL = 2kz — ky is a spectral parameter of the problem. 
System ([7]) is considered subject to the boundary conditions 

lim ao(x) = lim ai(x) = (8) 

corresponding to localized patterns in the x-direction. 

Now we establish several general properties of the solutions. First, multiplying (17a|) and 
( iTbl) respectively by dao/dx and dai/dx, adding conjugated equations, and integrating over 
X we obtain a necessary condition for existence of the localized modes: 

dV{x) 

, dx 

Next we notice that the Hamiltonian ([3]), corresponding to the effectively one-dimensional 
case we are interested in, can be written down in the form Huj = Hq + Hv{0) where 

Ho = j (^-^\ao,x\'^ - \ai^x\^ + doal + aodj^ dx (10) 

is the Hamiltonian of the two-wave interactions in the homogeneous medium and 

HviO = -\j V{x - {\a,\' + ^|ai|2^ dx (11) 

describes the effect of the periodicity. 

Assuming that a solution (ao(x), ai(x)) of Eq. ([7]) is given, we consider its infinitesimal 
shift C in the space, i.e. {ao{x — C),ai{x — Q). For the solution to be stable, the intro- 
duced shift must lead to increase of the energy. Thus for such a solution we obtain the 



(2|ao(x)p + |ai(x)p) rfx = 0. (9) 



conditions 



24| 



dHviO 



dC 



and 



dC 

C=o 



>0. (12) 

C=o 



The first one is nothing but Eq. ([9]) obtained above from the energetic arguments, while the 
second constrain acquires specific meaning for the definite symmetry solutions and will be 
considered below. 
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C. On phases of the stationary locahzed solutions 

In this section we prove that spatially localized solutions of Eqs. (j7[) are real. To this 

end, by analogy with 23|, l25|], we set aj{x) = pj{x)e'^^^^\ define 6 = 29i — Oq, and rewrite 
in the form 

Po,xx - Ol „,pQ - Vnpo + cos 6' = 0, (13a) 

+ ^o,xxPo + 2p?sin^ = 0, (13b) 

pi,xx - Ol^pi - ^Vnpi + 2pipo cos 6 = 0, (13c) 

26'i,^.pi,^. + Oi^^^pi - 2pipo sin 9 = 0, (13d) 

where Vq = V{x) — Q. Multiplying Eq. fll3dl) by pi and integrating we obtain 

d9i(x) 2 



PX 

/ p^(2;')po(x') sin 6'(x')(ix' (14) 
Jc 



p?(x) Jc 

where C is a constant. Substituting this last formula in Eq. fll3cp and considering the limit 
X ±oo together with the boundary conditions ([H]) we find out that a necessary condition 
for the existence of nonsingular solutions is that 



lim / pj{x')po{x')sme{x')dx' = 0. (15) 
Using the definition of pj and 6j this constrain can be rewritten as follows 

PX 

lim Im / dQ{x')a1{x')dx' 

x^±oo 

= J}fJ^£ ("^"("^'^^^-l^ + Vnix')\aoix')\'^ dx' 
= lim Re [ao(C)ao,2:(C) - ao(x)ao,a;(a:)] 

x^±oo 

= |po(C^)|%.(C) = 

where we have taken into account Eq. (fTa|) and performed integration by parts. From the 
last line of the above equalities we conclude that 6q does not depend on x, and thus without 
restriction of generality it can be chosen zero so that ao(x) can be chosen real. Then Eq. (I7ap 
implies that a\ is real. Thus one has two options: either oi is real or oi is pure imaginary. 
However, one readily concludes from the system (JTj) that the second option can be reduced 
to the first one by changing ^ —o-o- Therefore, in what follows, we will restrict our 
analysis to real stationary fields ao,i(a;). 
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D. Symmetry of solutions in even potentials 

For the sake of concreteness and to simplify the consideration below we set e{x) to have a 
definite parity. Then recalling the second of the constrains ( |T2l) we conclude that a solution 
which also has a definite parity, i.e. such that a'j{x) = a|(— a;), and strongly localized about 
X = 0, cannot satisfy (fT2l) if V{x) = —V{—x). Moreover, for strongly localized stable 
solutions Eg. (1121) is satisfied only if x = is a minimum of the periodic function V{x). Thus 
we can conjecture that stable solutions will be obtained centered about the minimum of the 
potential. 

Thus, we restrict further consideration to the case 



V[x) = V{—x), 



> (16) 

x=0 



dx'^ 

and consider a^ i either even or odd. 

As the next step we prove that among such solutions only even are allowed, or in other 
words we prove the following 

Proposition 1. If {aQ{x) , ai{x)} is a real solution of ^ with the potential [T^) . such that 
'^oi(^) ~ ^oii~^)' then necessarily: gq^x) = ao{—x) and ai{x) = ai{—x). 

Proof. First we prove that ao(x) is even. Assuming the opposite, i.e. that ao(a;) is an odd 
function and integrating (I Taj) with respect to x over the whole axis we obtains Pi = 0. 
Thus ai = and fl7a|l becomes linear and hence does not allow for the existence of spatially 
localized solutions. The contradiction we have arrived to, proves the claim that ao{x) must 
be even. 

Now we can prove that ai{x) is an even function. Again assuming that it is odd, what 

1 ci^^+^ai 

implies that ai(0) = 0, and designating (3n 



, , , ^ , -, 7^ the lowest nonzero 
(2n + l)! ^^Q^ 

derivative of ai(x) at x = 0, such that ai{x) = + o(x^""'"^) as x 0, from ( 17a|l 

we obtain that such a solution exists only if V{0) — Q = O(x^), and consequently from 
fl7a|) aQ^x) = 0{x^^"''^^^). Now one can estimate the orders of all terms in fl7bl) at x ^ 
as 0(x^"~^), 0(x^"+^), and 0{x^^'^^), i.e. fl7b|) cannot be satisfied in the vicinity of x = 0. 
We thus again arrive at the contradiction with the supposition about the existence of the 
solution. □ 
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E. On construction of spatially localized modes 

Let us now turn to the asymptotics of the solutions at |a;| oo. Since in this hmit, 
qj 0, it follows directly from (17b[) that ai{x) Ai{x) where Ai{x) is a decaying to zero 
solution of the Hill equation |22|] 



- ^ + i (^(^) - ^) ^1 = 0- (17) 

This equation has decaying (growing solutions) only if f2/4 belongs to a forbidden gap of 
the "potential" jV{x) (here we will not consider the cases where Q coincides with an edge 
of a gap). This last requirement we formulate as i7/4 G A^,!'' for some positive integer u'. 
Then it follows form the Floquet theorem that 

Ai(x)=Cie-'^^>i(a;) (18) 

where Ci is a constant, /ii is the corresponding Floquet exponent, defined by the frequency 
detuning to the gap, and 0i(x) is the real 27r-periodic function. 

Let us now turn to fITaj) . and notice that it can be considered as a linear equation for 
the unknown ao(a;). Taking into account, that for decaying solutions a{x), Q must belong 
to a gap 

A'-X of the potential V{x), i.e. n G a'-^] [ we emphasize that the frequency Q, 
here is the same as in Eqs. fl7b[) and f|T71) ]. we arrive at a necessary condition for existence 
of the localized modes: there must exist a nonzero intersection between gaps A^^>] and 



,(1) 

. (0) o ^ A (1) 



4:A^J for at least some integers z/' and u". Then, designating this intersection by A^, i.e. 



Ag = A^,; n 4:A^^/ and referring to it as a total gap, we conclude that if for a given f2 a 
localized mode exists, then Q E Ag. 
Now one can write down: 

ao(x) = ( A+(x) f A_{x')al{x')dx' + A_{x) j A+{x')a\{x')dx'\ . (19) 

where A± are the solutions of the eigenvalue problem 

+ {y{x) - ^]) A± = , (20) 



d^A. 



dx^ 

and W = A_^.{A_)x — {A^)xA_ is their Wronskian. Taking into account, that is in a 
gap of the spectrum of Eq. fl2U]) . we have that Aj-{x) = e^'^°^0±(x) with 0±(x) being real 
27r-periodic functions, and /iq being the respective Floquet exponent. Thus, in the case at 
hand W = -2/xo0+0- + 0+(0-)x - (0+)x0-- 
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Using the fact that J^^A+{x)al{x)dx = (this can be proved by substituting a1{x) 
from Eq. fl7a|) and integrating by part) the second term in (|T9l) can be rewritten in the same 
integration interval as the fist one 

ao{x) = ^ (^A+{x) A^{x')al{x')dx' - A_{x) A+{x')al{x')dx'^ . (21) 

F. Numerical study of localized modes 

Now one can use the shooting method to construct locahzed solutions of the system ([7j). 
As the first step one has to ensure that for a chosen structure the total gap exists. As 
the next step for a chosen VL G A^,, the Floquet exponents yUo,i and corresponding Bloch 
functions (pi and (j)± must be computed. Next, starting from some point x = xq far enough 
from the origin, where equation ( 17b[) is effectively linear, one can approximate the function 
ai(x) on the interval Xq < a; < oo by its linear asymptotics Ai{x)^ which is determined 
by (fTSi) with some small initial amplitude Ci. Substituting Ai{x) into equation (12T|) and 
computing the integrals one has to obtain the value of the function ao(xo) which now can be 
used to find ai {x) at the subsequent step of the spatial grid, x = xq — dx hy solving Eq. ( iTbl) 
(we use the Runge-Kutta method). Finally, by varying the shooting parameter Ci one has 
to satisfy the condition dao,i/dx = at the origin a; = (recall that we are looking for even 
solutions). 

Let us now show several possible solutions of the system ([7j). For the numerical simula- 
tions we assume the specific form of the dielectric permittivity: V{x) = — Vq cos(2x). In the 
left panel of FigiH an example of the band-gap structure of Eqs. (jlj) with a total gap (it is 
indicated by the shadowed region) is shown. In this particular case this is an overlap of the 
first gaps of each of the equation. The profiles of the corresponding solutions for different 
Q inside the total gap are shown in the right panels of Fig{Tl 

In all the pictures presented one observes that the amplitude of the signal is appreciably 
larger than the amplitude of the pump. Moreover, the amplitude of the mode itself (it 
is two-component in our case) decays and the pump amplitude become very small as the 
frequency approaches the bottom of the total gap. Notice that the top and bottom of the 
total gap coincide with the top of the gap for the pump signal and with the bottom of the 
gap for the signal. We also observe that the pattern of the pump is more sophisticated than 
the pattern of the signal component. 
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FIG. 1: In the left panel the band-gap spectra (q) for the pump (sohd Une) and 40]^ (q) for 
the signal (dashed line) are shown, v determines the corresponding band number. The shadowed 
domain shows the interval Ag = [1.89; 2.36] where the first gaps for the pump and the signal have 
overlapping and, consequently, where localized solutions of ^ can be found. In the right panels 
the profiles of the localized solutions ao{x) (solid lines) and ai(x) (dashed lines) with J7 € A^. Here 
Vo = 4. 

To investigate the stability of obtained solutions we integrate numerically the evolution 
equations ([T]) taking Ak = and starting with initial profiles presented in the right panels 
of Fig. [H In Fig. [2] we show the propagation of the localized modes for three different Q. 
We observe that the modes excited close to the bottom of the total gap are stable while the 
solutions in the vicinity of the top of the gap display unstable behavior. 
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FIG. 2: Dynamics of the modes, whose initial profiles are taken from the right panel of 
Fig. [1] and perturbed by 1% of their amplitudes. In [(a),(b)], [(c), (d)] and [(e), (f)] parameter 
n = 1.95; 2.19; 2.3, respectively. 

III. COHERENT STRUCTURES IN THE OPTICAL PARAMETRIC 

OSCILLATOR 

A. The model 



Now we proceed with the study of possible coherent patterns in an optical parametric 
oscillator and as the first step we recall derivation of the model equations, following 2l||. We 
consider the simplified case where in transverse direction there exists dependence only on one 
variable [we call it x: e = e{x)] and use the abbreviated notation 00,1(2^) ''") = cto.i(-2, x, r). To 
this end we impose the boundary conditions on ([1]), ([2]) assuming that the cavity is located 
on the interval < z < i: 



ao(0,r) = < + /2oe''°«o(£,r-T), 
ai(0,r) = Rie'^'ai{i,T -T) 



(22a) 
(22b) 



and allowing the amplitude of the driven field to depend on the transverse coordinates x, i.e. 
a™ = a^^{x). Here T = i/c is the delay time with c being the velocity of light, ^0,1 are the 
detuning parameters and -Ro,i are the reflectivity factors (in the following we will consider 
only the case of exact phase matching, Ak = 0). In the small amplitude limit, |ao,i| ^ 1, 
one can write down approximated solutions for pump and signal for the reduced Maxwell 
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equations ([T]) as the first order terms of the Mac-Laurin expansion 



ao{i,T) = ao{0,T) + -d^ao{0,T) + ieao{0,T) + iai{0,T) i, (23a) 
(0, r) + idlaiiO, r) + -eai{0, r) + 2zao(0, r)ai(0, r) £ . (23b) 
Assuming that the time variation of the solutions is very slow as compared to the delay 



T(9^ao,i(0) t) = "0,1(0^ r + T) - ao,i(0, r) 



(24) 



and considering limits of large reflectivity i?o,i ~ 1 and small detuning 6'o,i ^ 1 one comes 
to the model 



d^Ao = E{x) - (7 + tAo)Ao - -LqAo - Al 
drAi = -(1 + iAi)Ai - iLiAi + 2AoAi 



(25) 
(26) 



fsee also Refs. Refs. M and 



0, Refs. 3 



20], where the homogeneous counterpart of this system 



was derived respectively withing the framework of the multiple-scale expansion for excita- 
tions with smooth envelopes and as a continuum approximation for the cavity solitons in 
a periodic media, originally described by the coupled discrete equations). To derive these 
equations one uses the renormalized fields 

V{x) iaf 



i — iti 

rescaled time and space variables: 

T 



2 ' 1 - i?i 



(v4o,Ai,e,E), 



(27) 



(28) 



l-Ri ' 1-Ri 

as well as the constants 7 = {1 — Rq) / {1 — Ri) and Ao,i = 9o^i/ {Ri — 1). The linear operators 
£-0 1 have been defined in (111). 



B. Nonlasing states 

These states correspond to = and thus the respective Aq solves the stationary 
{dr = 0) equation 

Lo^'*^ + 2(Ao - 27)4'*^ + 2^^(a;) = 0. (29) 
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Let us consider the driven field having the same periodicity as the medium, i.e. E{x) = 
E[x + vr). Since the Bloch states constitute a complete orthonormal basis one can represent 

^(^) = E^-'?^S(^)- (30) 

(st) 

Looking for the respective stationary solution Aq also in a form of the expansion over the 
Bloch functions one readily obtains: 

4-"M=-2iE ,(.i^;;! . (31) 

u,q £V + 2(Ao - Z7) 

The obtained field pattern has an interesting particular limit which occurs when — 2Ao 
coincides with one of the band edges of the pump signal, i.e. when — 2Ao = £^±- Then 
A''q'^\x) = E(x)/'y, i.e. the pattern of the pump signal is the same as the pattern of the 
driven field (with accuracy of a factor I/7). 



C. Lasing states, particular case 1 



In the presence of a periodic potential of a general kind, finding explicit solutions for 
lasing states even in the simplest statement where Ai =const seems to be an unsolvable 
problem. A progress however is possible for periodic functions V{x) of a specific type. In 
order to show this we employ a kind of " inverse engineering" {23! and find an explicit form 
of V{x) assuming that Ai is a given complex constant: Ai = pe^'^ (here both p and (p are 
real constants). Then looking for the pump filed in the form Aq{x) = ao(a;)e^*'^, where ao{x) 
is a real function, we obtain from (l26l) 



ao[x) 



l+iAi + -V{x) 



which after substituting in fl25l) yields the system of two real equations: 

7 + 2, 



Ao + 7A1 + 



4 



-V{x) - 2Ei = 0, 



d^V{x) 
dx"^ 



- 2(Ao + 2Ai)\/(x) -r=^(x) 
+ 8(7 - AoAi + 2p2 - 2Er) 







(32) 



(33) 



(34) 



where Ej =Im(£'e ^*'^) and Er =Re(i?e ^*'^). Eq. viewed as an equation with respect 
to V{x), can be solved in quadratures: 



X 



^ 2(Ao + 2Ai)\/2 _ 16(7 - AoAi + 2p^ - 2Eu)V 



-1/2 



dV 



(35) 
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where without loss of generahty we have assumed that V{x) > with V^(0) = and \4(0) = 
0, i.e. that x = is a local minimum of the potential. This last condition implies existence 
of the threshold for Efj, which now must satisfy the inequality E^i > (7 — AqAi + 2p^)/2. 

The obtained expression ( l35l) depends on many physical parameters. Therefore, to sim- 
plify the analysis, we narrow the class of the potentials considering a case example of a 
one-parametric family of the periodic functions V{x): 

V{x) = Vo{k)sn'^ {f3{k)x, k) , (36) 

where sn(x, fc) is a Jacobi elliptic function with the modulus k, l3{k) = ^K{k), Vo{k) = 
6/c^/3^(fc) is the potential amplitude parametrized by k, and K{k) is a complete elliptic 
integral of the first kind. 

We observe that in practical terms the profile described by ( 136|1 is not too sophisticate: 
it can be reproduced with very high accuracy by only two or three harmonics for the elliptic 



parameter k not too close to one (see the discussion in |23| ) . 

Assuming that the refractive index of the cavity has the profile fl36l) . one finds that in 
order to support the stationary modes one has to apply the following form of the driven filed 

E{x) = e^'^ [Er + iEio + iEnsn^{l3x, k)] (37) 

where Ejq = (Aq + ■yAi)/2 and Eji = — (7 + 2)Vo(A;)/8 are fixed by the periodic structure 
and the cavity properties, while Eji is the control parameter which can be changed. Then 
the field pattern is obtained from (l32l) : 



p' = -lk'(3\k) - 1(7 - AoAi) + En, (38) 



«o(x) = ^ 



1 + iAi + ^Vosn^ {(3{k)x,k) 



(39) 



provided the elliptic modulus k is expressed in terms of the detunings Ao,i through the 
imphcit equation 

(1 + e)P\k) = -Ai - Ao/2. (40) 

The last equation means that Vq and Ao,i are not independent parameters and have to be 
matched and that the solutions we are dealing with exists only if at least one detuning is 
negative, or more precisely if 2Ai + Aq < 0. 
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Vo 



FIG. 3: The dependence Et^r '^•5 ^(^) for 7 = 1 and Ai = 1. 



Now, taking into account that p is a real constant one comes to the following threshold 
value Efhr for the control parameter 

Er > EtHr = lk^P\k) + 1(7 - AoAi). (41) 

As it is clear Ethr corresponds to p = 0. In Figj3]the dependence of the critical value of the 
driven field Et^r on the depth of modulation of the refractive index, Vq, is shown. 

In FigH] we show several patterns for different values of the amplitude of modulation of 
the refractive index, Vq, and of the control parameter, Er. Using initial dynamical equations 
(125|) . (!26|l we have checked that all the solutions presented are dynamically stable against 
initial small (of order of 5% of the intensity) perturbations. 



D. Lasing states, particular case 2 

The approach, based on the inverse engineering and described in Sec. IIII C[ can be gen- 
eralized. Indeed, let us look for a plane wave solution of the form 



^0 



e-''^\ Ai = Ai(x)e-*^" (42) 



where a; is a frequency to be determined and Ao^i{x) are real functions depending only on 
X (c.f. f l32|) where Aq = Ai/2 + V{x)/A). Here we use the form of the potential as in the 
Sec. II F, namely V{x) = — Vocos(2x). 
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FIG. 4: Patterns of the pump l^ol^ (solid line) and signal (dashed line) waves determined 

by (IMl), dSSD for different Vq and Er. In (a) Vq = 0.06, Er = 3, Aq = -4.03, in (b) Vb = 1.73, 
Er = 4, Ao = -4.88. in (c) Vq = 10.24, Er = 11, Aq = -9.63. The other parameters are 7 = 1 
and Ai = 1. 



Let us now require Aq to be a Bloch state of the linear eigenvalue problem 



-LqAq + AqAq = 2uAo . 



(43) 



As it is clear 2^01 must be a frequency bordering a gap, i.e. there must he 00 = jE,^^± + |Ao, 
since otherwise the eigenvalue Aq is complex. The amplitude of Ao{x) is a free parameter, 
so far, because (H3!) is linear. Let us now fix it by requiring the frequency u to border a gap 
of another linear eigenvalue problem, which reads 



^1 - 2Ao{x) 



Ai + AiAi = ujAi . 



(44) 



Then Ai is a Bloch state of this new problem where the effective periodic potential is given 
by V{x) = V{x)/A — 2Ao{x). Designating the gap edges of the Hill operator (l44l) by £[,2 
(notice that they depend on the amplitude of Aq{x)), we deduce that the amplitude of the 
pump signal Aq{x) is determined from the relation 



m + Ai 



lc(0) 



1 



Ao, 



(45) 



which can be solved numerically. The system fl25l) . fl26l) is satisfied if the complex field is 
chosen in the form 



E(x) 



- + Ai{x)+ti -Vix) - ^ + -Ao + ^Aoix] 



~2iujT 



(46) 
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FIG. 5: In (a) the band structure of the hnear eigenvalue problem (|44p and in (b) the profile of the 
corresponding effective potential, V{x), are shown. In [(c), (d)], [(e), (f)], and [(g), (h)] the solutions, 
Ai{x), corresponding to the first three edges of the bands denoted by points A (the edge Sq^), 
B (the edge Sq^), and C (the edge £^^) as well as corresponding profiles of the intensity of the 
driving field |i?(2;)|^ calculated from (j46p are shown. The profile of Aq{x) is calculated at the upper 
edge of the first band of the linear eigenvalue problem (|43p . The other parameters are 7 = 1 
and Vo = 1. 



In Figini we illustrate a band structure of (jUj) (panel a) as well as effective potential 
V{x) (panel b) with Aq{x) taken at the band edge and having the amplitude 1. The 
patterns of Ai{x) corresponding to the three lowest band edges denoted as A, B and C and 
the respective profiles of driving field |i?(x)p are shown in Fig. [5] (c)-(h). 



IV. CONCLUDING REMARKS 



In this paper we have carried out analytical and numerical studies of localized modes in 
optical parametric amplifiers and of periodic patterns in parametric oscillators generated by 
a properly chosen driving field in presence of diffraction and transverse periodic modulations 
of the refractive index. 

Using Bloch waves approach we have found that the presence of refractive index mod- 
ulations drastically affect the existence and properties of localized structures in the both 
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systems. In particular, spatial modulations lead to the formation of transverse phase locked 
(real fields) signal and pump in the form of a localized modes (gap solitons). Conditions 
required on the amplitude of modulations and external pump field for the existence of local- 
ized patterns in the parametric oscillators have been deduced. We have checked the stability 
of the obtained solution by means of integration of respective dynamical equations. It as to 
be emphasized however that the thorough analysis of the stability of gap solitons, which is 
complicated by the necessity of numerical generation of exact gap solitons, is left for further 
studies. In the present work we also did not discuss gap solitons of high order gaps, whose 
existence is strongly constrained (and may be even inhibited) by the necessity of the over- 
lapping of the spectra of two components, i.e. of the first and second harmonics (notice that 
the respective constrains do not exist in the theory of the single-component solitons in Kerr 
media). And meantime, the method elaborated for the particular case, where signal and 
idle waves were equal, allows straightforward generalization to the situation where resonant 
interaction among three different waves occurs. 
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